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c| , Recently a new general class of quark mass matrix ansatz has been 

proposed, which originates from some flavor symmetry. We extend 
that symmetry to the lepton sector and study the neutrino mass ma- 
trix and address the question of the baryon asymmetry of the universe 



X 

H 

C$ ' in this model. 



The question of flavor mixing and fermion masses can lead us to physics 
beyond the standard model, where several experimental results are present 
without any theoretical insight. As an attempt to derive relationship between 
the quark masses and flavor mixing hierarchies, mass matrix ansatz was 
suggested about two decades ago |J These ansatz hopefully will emerge from 
some definite theoretical consideration. 

Of the several ansatze of quark mass matrices, the canonical mass ma- 
trices of the Fritzsch-type have been generally taken to predict the entire 
Cabbibo-Kobayashi-Maskawa (CKM) matrix ||. However, this ansatz is 
ruled out because it predicts a top quark mass to be no longer than 100 GeV 
[§]. Recently this form of mass matrices has been generalized to accommo- 
date the large top quark mass while keeping the calculability of the ansatz 
||]. It was shown that this generalized mass matrix could originate from 
breaking of the maximal permutation symmetry 0, |£J. 

In this letter we study the consequence of the breaking of the maximal 
permutation symmetry in the lepton sector. Since this symmetry acts on both 
the left- and the right-handed quark fields, when we generalize it to the lepton 
sector, it would be natural to include the right-handed neutrinos. Then the 
same symmetry could act on the left- and right-handed leptonic fields in the 
same way as in the quark sector. We shall assume that lepton number is 
broken at some intermediate scale when the right handed neutrinos get a 
Majorana mass. Then the usual Higgs doublets will combine the left-handed 
neutrinos to the right-handed neutrinos through the Dirac mass term, so that 
the left-handed neutrinos get a small see-saw Majorana mass [0. There is 
no SU(2)l triplet Higgs scalar ||, which can break lepton number at some 
high scale and give a Majorana mass to the left-handed neutrinos. But, the 
maximal permutation symmetry will determine the form of the Majorana 
mass matrix for the right handed neutrinos and the Dirac neutrino mass 
matrix. From these matrices the see-saw mass matrix of the left-handed 
neutrinos will also be determined. With these neutrino mass matrices of 
the left-handed and the right-handed Majorana particles we shall study the 
possibility of generating a lepton asymmetry of the universe which will get 
converted to a baryon asymmetry of the universe before the electroweak 
phase transition [|, [H], [T]J . 

Let us begin by introducing the general form of the Fritzsch-type mass 
matrix in which a nonvanishing (2,2) element is introduced to accommodate 
a large top quark mass|4]]. This mass matrix form is achieved by breaking the 
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democratic flavor symmetry S(3)l'xS(3)r — > S(2)lxS(2)r — > S , (1)lxS , (1)/j. 
The resulting quark mass matrix takes the form, 




M q =\ADB\. (1) 



Although this matrix contains four independent parameters even in the case 
of real parameters, one can make additional ansatz to relate D to B in general 
so as to maintain the calculability ||. 

We shall now assume that the same flavor symmetry is also true in the 
lepton sector. In the S(3)l x S(3)r symmetric limit the charged lepton mass 
matrix has the form in the democratic basis: 




(2) 



which will lead to the largest contribution in the mass eigenstate basis. By 
making unitary transformation of the mass matrix Eq. (2) with the help of 



U — a/6 a/6 a/6 ' ( 3 ) 




one can obtain the diagonal mass matrix with only a nonvanishing (3,3) 
element m 33 = a. In order to account for the second generation masses, we 
break the S(3)l x S(3)r symmetry to S(2) L x S(2) R . This can be achieved 
by adding the following matrix to Eq.(2): 



(4) 



Although different choice of the parameters b and c is general, in this paper 
we set for simplicity, b = — c, which renders the mass matrix diagonal after 
unitary transformation with Eq.(3). We then get the second largest mass 
given by 777,22 = 26. Finally, the S(1)l x S(1)r symmetric mass matrix can 
be achieved by taking (1,2) and (2,1) elements to be nonzero (= d), which 




can then be rotated to the diagonal form with m\\ = d, 11122 = 2b + d and 
?"33 = a + c. 

We shall now work in the basis in which the charged lepton mass ma- 
trix is diagonal and real. For the S(3)l x S(3)r symmetry to be applicable 
we further consider that the model contains three right-handed neutrinos 
N a , a = 1, 2, 3 in addition to the usual left-handed neutrinos of three gener- 
ations. Furthermore, since the right-handed neutrinos are singlets under the 
standard model gauge group they can get Majorana masses at the tree level, 
whereas the left-handed neutrinos can not get any tree level Majorana mass 
since there are no SU(2)l triplet Higgs scalar. The left-handed neutrinos 
then remain light and get a see-saw mass @ due to their usual Dirac type 
coupling with the right-handed neutrinos. We can then write the neutrino 
mass matrix in the basis [un N aR ] as, 

M »=( °T I") ( 5 ) 

\rnjj M N J 

where each of these elements are 3x3 matrices. 

Assuming that the flavor symmetry under consideration comes from a 
transformation of the left- and the right-handed fields, the right-handed Ma- 
jorana mass matrix will have a S(3)r x S(3)r symmetry to start with, which 
will then break to S(2) R x S(2)r and then to S(1)r x S(1)r in succession. 



The S(3)r x S(3)r symmetry gives a mass matrix of the form [12 



A B B\ 
M N =\B A B (6) 

B B A) 

in which we can choose B = without loss of generality. Otherwise, we 
can keep the B and diagonalize it with the unitary matrix with the first row 
and the third row interchanged from Eq.(2) to a form proportional to a unit 
matrix M^ — A I. Then the combination Ni + N 2 + N 3 gets a mass, which we 
identify with the largest mass of the right-handed neutrino iV e ij. Identifica- 
tion of the right-handed electron neutrino as the heaviest one is neccessary to 
explain the large mixing between the left-handed \x and r neutrinos, which is 
required to explain the atmospheric neutrino anomaly. We shall discuss this 
point in details at some later stage. In the same way, the S(2)r x S(2)r sym- 
metry will contribute to the (2,2) diagonal elements by an amount C, and the 



Table 1: Present experimental constraints on neutrino masses and mixing 



Solar Neutrino [|l3j 
(Small angle MSW) 


Am 2 ~ (0.5 - 1) x l(T 5 e1/ 2 
sin 2 20 ~ KT 2 - 10~ 3 


Atmospheric Neutrino 14 


Am 2 ~ (0.5 - 6) x 10~ 3 e1/ 2 
sin 2 26 > 0.82 


Neutrinoless Double Beta Decay [[L5f 


m Ve < QAQeV 


CHOOZ [|1§ 


Am 2 eX < 10- 3 eV" 2 
(or sin 2 26 eX < 0.2) 



corresponding state may be identified as the N^r. Finally, the S(1)r x S(1)r 
symmetry will contribute only to the (3,3) element by an amount D. The 
resulting matrix in the basis [N e R N^r N t r] then becomes, 



M, 



N 



diag[A, C, D]. 



(7) 



We now assume that as in the quark sector the largest contribution to the left- 
handed neutrino mass comes in the 5(3)^ x S(3)r symmetric limit, then the 
heavy right-handed neutrinos follow an inverted hierarchical pattern M^ e ^> 
M^fj, 3> M^ T and N t r becomes the lightest right-handed Majorana neutrino. 
As we shall demonstrate, this inverted hierarchical pattern can explain the 
present experimental results of the neutrino masses (which is given in table 

As it is well known, it is not possible to explain the atmospheric neutrino 
anomaly [fPH , solar neutrino problem |13| and the LSND result |17| simul- 
taneously in a three generation model, since the mass squared differences 
required in the three cases are widely different. One needs at least one more 
light sterile neutrino state which mixes with the ordinary neutrinos. On the 
other hand, such sterile neutrinos are severely constrained by the present 
limit on the nucleosynthesis bounds [|llj, which can hardly accommodate a 
fourth neutrino. So, in our analysis we consider only three light neutrinos 



and thus do not try to accommodate the LSND result. There are three al- 
lowed regions of the parameter space for explaining the solar neutrino data 



8], out of which we consider only the small angle MSW solution, which 
comes out naturally in this model. 

The structure of the Dirac mass matrix in the neutrino sector will be 
similar to that of the quark sector and is given by, 

(8) 




However, the hierarchy in the different elements could, in general, be dif- 
ferent. So, we try to determine the different parameters from experimental 
inputs, rather than justifying them from some theoretical reasonings. 

At this point we can justify the requirement for the inverted hierarchy of 
the right-handed Majorana neutrino masses. Thanks to the largely hierarchi- 
cal right-handed Majorana neutrino mass matrix, the main contribution to 
the left-handed Majorana masses comes from the lightest right-handed neu- 
trinos in the see-saw mechanism. If N e R becomes the lightest right-handed 
Majorana neutrino, then the largest element of the effective left-handed Ma- 
jorana mass will be the [/i/x] element only. On the other hand, the atmo- 
spheric neutrino anomaly requires a large mixing between v^ and u T , which 
requires the [/it] and [rft] elements to be comparable or larger than the [/i/i] 
element. Thus this cannot explain the large v^ and v T mixing. If we assume 
that N t r is the lightest right-handed Majorana neutrino and y and z are of 
the same order of magnitude pOfl , then all the four elements [/i/x], [/xr], [r/i] 



and [tt] will be of the same order of magnitude, which ensures maximal 
mixing between v^ and v T ]2T| . 

Along with the inverted hierarchy condition, M/y e ^> M^ 3> M^ T , we 
also assume that the elements of the Dirac mass matrix are not largely hi- 
erarchical. Then the largest contribution to the effective 3x3 mass matrix 
of the left-handed neutrinos will come from the see-saw contribution from 
Mtv t and then the next leading contribution will come from M^ . Then the 
effective 3x3 neutrino mass matrix can now be written as 

m v = m T D M^ l m D 



'2 



/ x z xt xy \ 

/ M N/i M Nfi M N)i \ 

xt y 2 yz 

M N ^ M Nt M Nt 

xy yz z 2 . 

\ Mjv m M Nt M Nt ) 



(9) 



where we have ignored the next to leading terms of order 0(1/ M^) and 
0(1/Mff ). This means that the mass squared difference between v^ and v T 
vanishes. So, for the nonvanishing mass squared difference we need to keep 
the higher order terms. However, to calculate the mixing matrix, we can 
ignore the next to leading terms without loss of generality. 

Since we are interested in the maximal mixing solution for v^ and v T and 
the small mixing solution for v e and u^, the neutrino mixing matrix can be 
parameterized by 



U v 



1 N 


\ I 


C\ 


Si 





c 2 s 2 


' 


-si 


Cl 





s 2 c 2 , 


' V 








f 


Ci s 1 


o \ 


i 






Sl Cl 

si Cl 
V2 \/2 


71 


) 







(10) 



where ci = cos$i,C2 = sin 9 2 and we have taken c 2 = s 2 = l/\/2. Then, 
the favored solution for the solar and atmospheric neutrino anomaly leads 
to ci ~ 1, si ~ 0.05, m 2 ~ 3 x 10~ 3 eV and m 3 ~ 3 x 10~ 2 eV. From these 
results, one can approximately estimate the neutrino mass matrix 



m v = 



U u ■ diag[m 1 , -m 2) m 3 ] ■ XP V (11) 

' 10 -6 10 -4 1Q -4 \ 

10~ 4 1.5 x 10~ 2 1.5 x 10~ 2 (12) 

v 10~ 4 1.5xl0~ 2 1.5 xlO" 2 / 

In this case, by choosing y and z to be of the same order of magnitude 
we can ensure the maximal mixing between v^ and v T . To be precise, the 
large mixing between v^ and is T , as observed in the atmospheric neutrinos 
at SuperKamiokande(i.e., sin 2 29 > 0.82), implies 0.64 < y/z < 1.56. On 
the other hand, the ratio x/y ~ x/t gives the mixing between v e and v^ and 
could be very small. For a reasonable choice, x/y ~ x/t ~ 10~ 2 , the solar 



neutrino anomaly can be explained by the small angle MSW solution, which 
is one of the favored solutions for the solar neutrino anomaly 



In all the models for neutrino masses, another related question remains to 
be answered is the generation of the baryon asymmetry of the universe. The 
Majorana mass term of the right handed neutrinos introduces lepton number 
violation and hence (B — L) violation. If this interaction is slower than the 
expansion rate of the universe when the right handed neutrinos decay (at 
T = Mjv) and there is enough CP violation, this interaction can generate a 
lepton asymmetry of the universe. At finite temperature above the critical 
temperature of the electroweak phase transition sphaleron processes are in 
thermal equilibrium and (B+L) number violating interactions are very fast 
HJ . This will then relate the lepton asymmetry or the (B — L) asymmetry 



generated during the right handed neutrino decay to the baryon asymmetry 
of the universe before the electroweak phase transition. This remains to be 
the most interesting scenario for the understanding of the baryon number of 
the universe, which is referred to as leptogenesis fllOl , |22| . 

Since lepton number violation is the source of leptogenesis, it also de- 
pends on models of neutrino masses. It has recently been argued |23| that 



the see-saw mechanism of neutrino masses is the most preferred one for the 
generation of the baryon asymmetry of the universe in supersymmetric infla- 
tionary models. We shall now see if this could be implemented in the present 
scenario. 

In the see-saw mechanism for neutrino masses the right-handed neutrino 
decay could generate a lepton asymmetry of the universe, which then gets 
converted to a baryon asymmetry of the universe. CP violation comes from 
an interference of the tree level diagrams with the self-energy type and vertex 



correction type diagrams [jTy, p2|. We start with the lagrangian 



£ = MiNc R N iR + h ai £ aL <pN lR + h.c, (13) 

where £ a L are the light leptons, <p is the usual Higgs doublet of the standard 
model, h ai are the complex Yukawa couplings and a is the generation index. 
We have chosen the Majorana mass matrix of the right-handed neutrino 
to be real and diagonal. Then the decay of the lightest neutrino N t r can 
generate a lepton asymmetry which can then generate a baryon asymmetry 
of the universe before the electroweak phase transition. The decay of the 
heavier neutrinos can also generate an asymmetry, but that asymmetry will 
be washed out before the lightest right-handed neutrino decay. 
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Due to the Majorana masses of the right-handed neutrinos, their decay 
violates lepton number, 

N iR -> 4i + f (14) 

+ 0. (15) 



9 c 



CP violation comes from an interference of the tree level diagram for the 
decay of Nm with the one loop diagrams shown in figure 1. There are two 
contributions coming from the interference of the tree diagram with the one 
loop vertex correction and another one with the one loop self energy dia- 
gram. When the masses of the heavy neutrinos are degenerate, there can be 
large contributions coming from an interference of the self energy diagrams. 
However, in the present case when the masses of the right handed neutrinos 
are hierarchical, the contribution coming from the two sources are equal. In 
this case the contributions from both the diagrams add up to give the total 
amount of CP violation. First, the heaviest right handed neutrino will decay 
and then the second heaviest, when it may or may not generate any lepton 
asymmetry. When the lightest one (N tR ) decays, it will first erase any pre- 
existing asymmetry and then generate the final lepton asymmetry, which is 
presented by 



£a Wfl - W) " £« Wfl 



&* 



E a r(N TR -> 4L0t) + £ a t(n tR - t aL <t>) 
3 imfEJ/CM EMM)} m 2 Nt 

where i — e, fj, and I(x) = x 1 / 2 [l + (1 + x) ln(x/(l + x))\. 

Similar to the Jarlskog invariant for CP violation in the quark sector, in 
this case a different combination enters, which has its origin in the Majorana 
nature of the neutrino masses. The Majorana nature of the right handed neu- 
trinos will imply that there are new Majorana phases which will contribute to 
the CP violation. In the basis we are working, where the right handed mass 
matrix is real and diagonal, all these Majorana phases has been transferred 
to the elements of mo and the quantity J2 a (ha T h a i) becomes a rephasing 
invariant quantity. In the limit x — 0, this becomes equivalent to a two 
heavy neutrino scenario. In this case after considering the overall rephasing 



there remains only one CP phase, which vanishes when t = 0. Thus, in the 
Fritzsch type of mass matrix, the CP violation for the generation of the a 
baryon asymmetry of the universe will be highly suppressed. 

In terms of the elements of trd, the amount of CP violation is given by, 



(t 



16tt 



lm(y*t + z*yf ( ' M N \ Im(y*x) 2 ( M, 



(M 2 + M 2 ) \m n J (M 2 + M 2 ) \m. 



l N T 



L N e 



(17) 



We can now choose the overall phase so that z*y is real. Then for t = 0, 

the first term vanishes and assuming that the CP phases are all similar 

t A/f 
the amount of CP violation gets suppressed by an amount, < 10~ 4 . 

yM Nr 

With this additional suppression it will be impossible to explain the baryon 
asymmetry of the universe. Thus, for the generation of the lepton asymmetry 
of the universe, it is extremely important that we make the [2,2] element 
non- vanishing, which required to accommodate the heavy top quark mass in 
Fritzsch type (mass matrix) ansatz. 

For the generation of a lepton asymmetry of the universe one more in- 
gredient is neccessary, namely, the out-of-equilibrium condition. Whether a 
system is in equilibrium or not can be understood by solving the Boltzmann 
equations. But a crude way to put the out-of-equilibrium condition is to say 
that the universe expands faster than the interaction rate J23J]. This may be 
stated as 

T N <H=1.7Jg;— (18) 

where Tn is the interaction rate under discussion, g* ~ 10 2 is the effective 
number of degrees of freedom available at that temperature T, and Mp is 
the Planck scale. 

In the case of right-handed neutrino decay, the asymmetry is generated 
when the lightest right handed neutrino decays. Before its decay, the pre- 
existing lepton asymmetry, if any, is washed out by its lepton number vio- 
lating interactions. Just before the lightest right handed neutrino decays, it 
satisfies the out-of-equilibrium condition 



\h 



ol 



12 



-M N < H at T = M N (19) 

16tt V ; 

if the masses of the heavy neutrinos are larger than M^ > 10 7 GeV. We 
consider that the reheating temperature is not too high so that after reheating 
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gravitinos are not produced in large numbers which can then overdose the 
universe [25|. This implies that for the lightest right-handed neutrinos to 
be produced after reheating, we must have M Nr < 10 10 GeV. This reduces 
the uncertainty in the scale of the heavy neutrinos to some extent. The 
explanation of the atmospheric neutrino anomaly, i.e., the mass of v T to be 
around 10~ 2 then requires, y, z, t ~ 0.1 GeV for M Nt ~ 10 8 GeV. Taking 
the vacuum expectation value of the Higgs doublet field to be around 100 
GeV, we get a lepton asymmetry of the universe to be around, til ~ — ~ 
10 -8 sin 5, where 5 is the CP phase in the couplings of the Dirac neutrino 
mass matrix rnp. So, an CP violating phase of the order of sin 5 ~ 0.01 
can generate enough lepton asymmetry of the universe. During the period 
10 12 GeV > T > 10 2 GeV the sphaleron transitions will be very fast and this 
lepton asymmetry will get converted to a baryon asymmetry of the universe, 

In summary, we extended a new general class of quark mass matrix ansatz 
to the leptonic sector to obtain the neutrino mass matrix. We showed that an 
interesting inverted hierarchical pattern for the heavy right-handed neutrinos 
can accommodate the atmospheric neutrino oscillation with the maximal 
v^ — * v T mixing and the small angle MSW solution of the solar neutrino 
deficit. It turned out that the baryon asymmetry of the universe comes out 
to be correct for this particular form of neutrino mass matrices. 

Acknowledgments One of us (US) would like to thank the organisers of 
the journal club and the Korea Institute for Advanced Study for hospitality, 
where this work was initiated. 



10 



References 

[1] S. Weinberg, in a Festschrift for I. I. Rabi, ed. L. Motz (N.Y. Academy 
of Sciences, N.Y. 1977); A. deRejula, H. Georgi, and S. L. Glashow, 
Ann. Phys. (N.Y.) 109, 258 (1977); F. Wilczek and A. Zee, Phys. Lett. 
B79, 418 (1977); H. Georgi and D. V. Nanopoulos, Nucl. Phys. B155, 

52 (1979). 

[2] A. C. Rothman and K. Kang, Phys. Rev. Lett. 43, 1548 (1979); Phys. 
Rev. D24, 167 (1981); H. Fritzsch, Phys. Lett. B73, 317 (1978); Nucl. 
Phys. B155, 189 (1979). 

[3] K. Kang and S. Hadjitheodoridis, Phys. Lett. B193, 504 (1987); K. 
Kang, J. Flanz, and E. Paschos, Z. Phys. C 55, 75 (1992) See also H. 
Harari and and Y. Nir, Phys. Lett. B 195, 586 (1987). 



[4] K. Kang and S. K. Kang, Phys. Rev. D 56, 1511 (1997); [Kep= 
ph/980233C| . 



[5] H. Fritzsch and Z. z. Xing, Phys. Lett. B 413, 396 (1997); [hep 



ph/9904286| ; A. Mondragon, E. Rodriguez- Jauregui, Phys. Rev. D 59, 
093009 (1999); |hep-ph/9906429| ; D. Falcone, frep-ph/99053161 M. Rand- 



hawa, V. Bhatnagar, PS. Gill, M. Gupta, |hep-ph/9903428] 



[6] See also, S. L. Adler, Phys. Rev. D 59, 015012 (1999); [hep-ph/9711393 



P. S. Gill and M. Gupta, Phys. Rev. D 56, 3143 (1997); K. Kang, S. K. 
Kang, C. S. Kim and S. M. Kim, |hep-ph/98084T9| . 

[7] M. Gell-Mann, P. Ramond and R. Slansky, in Supergravity, eds. P. 
van Nieuwenhuizen and D. Freedman, (North-Holland, 1979) p. 315; T. 
Yanagida, in Proc. Workshop on Unified Theories and Baryon Number 
in the Universe, eds. A. Sawada and A. Sugamoto (KEK, 1979), p. 95; 
R.N. Mohapatra and G. Senjanovic, Phys. Rev. Lett. 44, 912 (1980). 

[8] E. Ma and U. Sarkar, Phys. Rev. Lett. 80 (1998) 5716; U. Sarkar, 
Phys. Rev. D 59 (1999) 031301; E. Ma. M. Raidal and U. Sarkar, [hei> 
ph/98ll24T]| ; Eur. Phys. J. C8 (1999) 301. 

[9] E. W. Kolb and M. S. Turner, 1989, The Early Universe (Addison- 
Wesley, Reading, MA). 

11 



[10] M. Fukugita and T. Yanagida, Phys. Lett. B 174, 45 (1986). 

[11] S. Yu. Khlebnikov and M.E. Shaposhnikov, Nucl. Phys. B 308, 885 
(1988); J.A. Harvey and M.S. Turner, Phys. Rev. D 42, 3344 (1990). 

[12] K. Kang, S. K. Kang, J. E. Kim and P. Ko, Mod. Phys. Lett. A 12, 
553 (1997); K. Kang and S. K. Kang, |hep-ph/9802328| ; S. K. Kang and 
C. S. Kim, Phys. Rev. D 59, 091302 (1999); H. Fritzsch and Z. z. Xing, 
Phys. Lett. B 372, 265 (1996); [hep-ph/ 9807234] ; Phys. Lett. B 440, 



313 (1998); ; M. Tanimoto, |Kep-ph/9807517l ; M. Fukugita, M. Tanimoto 
and T. Yanagida, Phys. Rev. D 57, 4429 (1998). 

[13] Super-Kamiokande Collaboration : Y. Fukuda et al, Phys. Rev. Lett. 
81, 1158 (1998); Talk by Y. Suzuki at Neutrino'98, Takayama, Japan 
(1998). 

[14] Super-Kamiokande Collaboration : Y. Fukuda et al, Phys. Rev. Lett. 
81 (1998)1562 ; |hep-ex/98050"06| ; Phys. Lett. B433 (1998) 9; T. Kajita, 



|hcp-cx/981000T 



[15] Heidelberg-Moscow Collaboration : H.V. Klapdor-Kleingrothaus, in 
Proc Lepton and Baryon number violatoin, Trento, April 1998; M. 
Giinther et al, Phys. Rev. D 55 (1997) 54; L. Baudis et al, Phys. Lett. 
B 407 (1997) 219. 

[16] CHOOZ Collab., M. Apollonio et. al, |hep-ex/971100l . 

[17] C. Athanassopoulos et. al, Phys. Rev. Lett. 75 (1995) 2650. 

[18] N. Hata and P. Langacker, Phys. Rev. D56 (1997) 6107; J.N. Bahcall, 
P.J. Krastev and A.Yu. Smirnov, |hep-ph/9807216j . 

[19] V.A. Kuzmin, V.A. Rubakov and M.E. Shaposhnikov, Phys. Lett. B 
155 (1985) 36. 

[20] E. Ma, D.P Roy and U. Sarkar, Phys. Lett. B444 (1998) 391; E. Ma 
and D.P. Roy, Phys. Rev. D 59, (1999) 097702. 

[21] R. Barbieri, L.J. Hall and A. Strumia, |hep-ph/980833l . 



12 



[22] P. Langacker, R. Peccei and T. Yanagida, Mod. Phys. Lett. A 1, 541 
(1986); A. Acker, H. Kikuchi, E. Ma and U. Sarkar, Phys. Rev. D 48, 
5006 (1993); M. Flanz, E.A. Paschos and U. Sarkar, Phys. Lett. B 345, 
248 (1995); L. Covi, E. Roulet and F. Vissani, Phys. Lett. B 384, 169 
(1996); M. Flanz, E.A. Paschos, U. Sarkar and J. Weiss, Phys. Lett. 
B 389, 693 (1996); W. Buchmuller and M. Plumacher, Phys. Lett. B 
389, 73 (1996); Phys. Lett. 431, 354 (1998); J. Faridani, S. Lola, P.J. 
O'Donnell and U. Sarkar, Eur. Phys. J. C 7, 543 (1999); A. Pilaftsis, 
Phys. Rev. D 56, 5431 (1997); Nucl. Phys. B 504 (1997) 61. 



D. Delepine and U. Sarkar, |hep-ph/981iTT9l ; E. Ma, M. Raidal and U. 
Sarkar, |hep-ph/990i406 . 



[24] J.N. Fry, K.A. Olive and M.S. Turner, Phys. Rev. Lett. 45 (1980) 2074; 
Phys. Rev. D 22 (1980) 2953; Phys. Rev. D 22 (1980) 2977; E.W. Kolb 
and S. Wolfram, Nucl. Phys. B 172 (1980) 224. 

[25] J. Ellis, G. B. Gelmini, J. L. Lopez, D. V. Nanopoulos and S. Sarkar, 
Nucl. Phys. B373 (1992) 399; M. Boltz, W. Buchmuller and M. 
Plumacher, Phys. Lett. B443 (1998) 209; D. J. H. Chung, E. W. Kolb 
and A. Riotto, [hep-ph/ 98094531 . 



13 



N 




N ,\ N 




4> 



N / 
/ 



N 



i L c 




\ 
4> 



Figure 1: Tree level and one loop vertex and self energy diagrams for the 
generation of lepton asymmetry of the universe 
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